The head leads the body: a curvature-based kinematic description of C. elegans 
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Caenorhabditis elegans, a free-living soil nematode, propels itself by producing undulatory body 
motion and displays a rich variety of body shapes and trajectories during its locomotion in complex 
environments. Here we show that the complex shapes and trajectories of C elegans have a simple 
analytical description in curvature representation. Our model is based on the assumption that the 
curvature wave is generated in the head segment of the worm body and propagates backwards. We 
have found that a simple harmonic function for the curvature can capture multiple worm shapes 
during the undulatory movement. The worm body trajectories can be well represented in terms of 
piecewise sinusoidal curvature with abrupt changes in amplitude, wavevector, and phase. 



AUTHOR SUMMARY 

A thorough understanding of the kinematics of C. el- 
egans is crucial to investigate complex behaviors such as 
locomotion, chemotaxis and thermotaxis. In this work, 
we show that C. elegans propels itself and turns using a 
simple set of elementary movements, each of which cor- 
responds to a single harmonic curvature wave. The cur- 
vature wave that originates in the head segment propa- 
gates along the worm body, and abrupt changes in the 
amplitude, wavevector, and phase of the wave result in 
the corresponding changes of the overall direction of a 
moving nematode. We demonstrate that our model can 
be used not only to describe the individual worm body 
postures but also to determine the entire trail with its 
multiple shapes. A deeper insight into the repertoire of 
simple elementary movements of C. elegans will advance 
qualitative studies of nematode behavior. We also ex- 
pect that it will enable us to decipher the neural-sensory 
control of the nematode. 



INTRODUCTION 
Motivation 

A one-millimeter long soil dwelling nematode 
Caenorhabditis elegans is a model organism for compre- 
hensive investigations spanning genetics [l|, neural 
control [sl, sensory transduction [5[, behavior [gI, 
and locomotion With a total of 959 nongonadal 

cells, of which exactly 302 are neurons [3], C. elegans is 
a very simple organism. Yet, it can efficiently move in 



complex environments [1 
behavior to mechanical [14 



purposefully adjust its 
, chemical [17., ^], and 



thermal stimuH [Il-il, and learn 

So far many studies of nematode behavior have focused 
on gross characteristics, such as the probability of turns 
or the fraction of worms that approach a chemical signal 
within a prescribed time [2J-[26|. Also, typical body pos- 



tures of C. elegans (e.g., sinusoidal undulations, shallow 
turns, and sharp Q- or loop-turns) have been classified 
and used to characterize worm trajectories. Such descrip- 
tions of worm motion, however, are inadequate because 
they do not address more advanced questions regarding 
worm behavior: for example, how the nervous, sensory, 
and motor systems of C. elegans control its body to gen- 
erate efficient propulsion and move the nematode towards 
a food source or away from danger. To gain insight into 
such complex behaviors a quantitative approach to model 
worm dynamics is a prerequisite. 

An important contribution to develop quantitative 
models of worm gait has been made by Stephens et al. 
[27]. They have shown that a set of complex shapes that 
C. elegans assumes during crawling has a simple rep- 
resentation in orientational coordinates intrinsic to the 
worm body. In their model a significant fraction of such 
shapes are expressed in terms of several eigenmodes of a 
correlation matrix for orientation of body segments. This 




FIG. 1. Trail and body postures of C. elegans crawling 
on agar. The figure shows two superimposed snapshots of 
the nematode crawling on a 4 wt% agar surface. 
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FIG. 2. Schematic of a C. elegans trail. The worm body 
is a portion of the entire trajectory. 



" eigenworm" description provides a low dimensional rep- 
resentation of the complex space of worm postures, which 
is a crucial step towards developing realistic models of the 
worm behavior. 

On the other hand, the eigenworm picture [27] is not 
intuitive, and it misses an important characteristic of 
worm trajectories, recently described by Kim et al. [28^ . 
Based on an analysis of trajectory images, they have pro- 
posed [28] that shallow turns of crawling worms can be 
described using a model where the worm performs simple 
sinusoidal undulations with suddenly changing parame- 
ters such as wavelength and amplitude. While this model 
agrees well with experimental observations of C. elegans 
performing low-amplitude undulations [28], it is not suf- 
ficient for a description of 1^- turns and other sharp turns. 



Piecewise-harmonic-curvature model 

Our curvature-based description of C. elegans motion 
is motivated by geometrical observations of trails of ne- 
matodes moving without slip on a soft agar substrate. In 
this work we build on the ideas of [27] and [28] to develop 
a simple curvature-based description of worm kinemat- 
ics, obtaining a model that overcomes the shortcomings 
described above. As in [28|, we analyze entire trails of 
crawling nematodes (rather than their individual body 
shapes) to capture sudden changes of the gait that the 
worm uses to navigate its environment. However, we 
formulate our description in terms of intrinsic geometric 
quantities (arc length s and local curvature n)^ rather 
than real-space coordinates, and this allows us to go be- 
yond describing only shallow turns and low-amplitude 
undulations. As a result, we obtain a low-dimensional 
representation of all worm postures and trajectories in 
terms of uncomplicated basis functions. 

A typical trail is shown in Fig. (TJ where two images 
are superimposed to depict the worm postures at differ- 
ent times. The worm makes deep U-shaped undulations 
throughout the trajectory. The depth of the undulations 



varies in time, and due to these variations the overall 
direction of the worm motion changes. The worm occa- 
sionally assumes an (^-like shape (as seen in the top right 
corner of Fig. [1]), which results in a significant change of 
direction. 

Assumptions. Our piecewise-harmonic-curvature 
(PHC) model is built on the following three assumptions: 
(i) the worm controls its body posture by propagating 
the curvature backwards along the body, (ii) the curva- 
ture wave is generated in the front (head) segment of 
the worm body; the time variation of the curvature of 
this segment serves as the boundary condition for the 
curvature wave and determines this wave completely, 
and (iii) the variation of the head segment curvature can 
be described in terms of a simple piecewise-harmonic 
function. 

As shown below, our model is confirmed by a detailed 
analysis of trails of worms crawling without slip on agar 
substrate. Trail segments described by harmonic varia- 
tion of the curvature are clearly identified, and the jumps 
of the curvature- wave amplitude, wave vector, and phase 
from segment to segment are determined. We thus ob- 
tain an analytic representation of the entire trail in terms 
of a relatively small number of parameters. 



RESULTS 

Curvature representation of the worm trail 

It is convenient to describe nematode trails and instan- 
taneous body shapes using the parametric description 



x = x{s), y = y{s) 



(1) 



for the Cartesian x and y coordinates of all points along 
the trail. The parameter s in Eq. ([1]) represents the arc 
length along the trail. The worm body posture at time t 
corresponds to a portion of the trail between the points 



St = ST{t), 5h = ST{t) + L, 



(2) 



where st is the tail position, su is the head position, and 
L is the length of the worm. If the worm is moving with 
a constant velocity Vg^ the tail position is given by 



ST{t) =ST{0)^Vst. 



(3) 



The above description of the trail geometry is illustrated 
in Fig. El 

The shape of the worm trail in the parametric form ([1]) 
can be evaluated from the local curvature = hz{s) using 
the Frenet-Serret equations 



di{s) 

ds 



K.{s)h{s), 



where 



dx ^ dy ^ 



dx ^ dy ^ 



(4) 



(5) 
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FIG. 3. Family of shapes obtained from a harmonic 
curvature function. Solid lines show the solutions of 
Frenet-Serret equations JB)) with the harmonic curvature (|8]) 
with different values of normalized amplitude Ajq. The 
dashed lines correspond to the real-space sinusoidal form Q. 



are the unit vectors tangent and normal to the trail, re- 
spectively. Combining relations dH and (|5|) yields a set 
of coupled second-order differential equations 



d\ 
d^ 
d^ 

d52 



d^ 

d5 

dx 
d^ 



(6a) 
(6b) 



for the trail shape ([T]) . Equations (j6j) are solved with the 
initial conditions 



x(^s = 0) = xo, 
y(s = 0) = ?/o, 



dx(5) 



As 
d^(5) 



s=0 



ds 



(7a) 
(7b) 



s=0 



where (xo,yo) is the position of the worm tail at the 
beginning of the trail and the tangent unit vector, to = 
i'^xo^ty^) describes the orientation of the tail at t = 0. 



FIG. 4. Body postures of C. elegans crawling on agar. 

Left column: Experimental images. Right column: corre- 
sponding fits of the harmonic-curvature representation ([8]) 
(solid lines) to the skeletonized data (open circles). For clar- 
ity, only every fifth point of the skeleton data is plotted in all 
the figures. Figure (B) also shows the real-space sinusoidal fit 
dS]) (dashed line), which significantly deviates from the skele- 
ton data. 



Harmonic-curvature variation 

We argue that worm trails (and the corresponding se- 
quences of worm shapes) can be accurately described in 
terms of a simple harmonic-curvature variation 



hz{s) = Asm{qs - 



(8) 



where A is the amplitude, q is the wave vector, and (j) is 
the phase of the harmonic wave. Below we demonstrate 
that a combination of segments of the lines obtained by 
integration of the curvature (|8|) with different values of 
the parameters A, and (j) can be used to construct the 
entire worm trail. 

The family of shapes obtained by changing the ampli- 
tude A in Eq. ([8j), with q and (j) fixed, is depicted in Fig. 
[3l Other shapes, corresponding to different values of q 
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FIG. 5. Body posture of C. elegans with two distinct curvature modes. (A) The experimental image of the worm. 
(B) The skeleton data (open circles) and the harmonic-curvature fit for the tail segment (solid line). The extension of the fit 
that does not follow the head segment is shown by dashed line, and the point where the curvature mode changes is indicated 
by a filled circle. (C) The best fit of the two- mode harmonic-curvature representation (|8]) (solid line). Since the line is obtained 
by integrating second-order differential equations (|6|), it is continuous and has a continuous slope. The insets in (B) and (C) 
show the local fit error (|17p along the skeleton of the nematode. For the single-mode fit (B) the error rapidly increases after the 
point indicated by the filled circle, whereas for the continuous two-mode fit (C), the local error is below 1% along the whole 
body. 




FIG. 6. Loop turn with two distinct curvature modes. 

(A) Experimental image. (B) The two-mode harmonic- 
curvature representation (solid line); the skeleton data (open 
circles); the location of the mode change (filled circle). 

and can be obtained from this family by affine trans- 
formations. A comparison of the lines plotted in Fig. [3] 
with the worm tracks and body shapes shown in Figs. [1] 
and [4]-[9] reveals close similarity. For example, the shape 
of the line with the dimensionless amplitude A/q = 2 
closely resembles the 1^- turns made by the worm in Figs. 
[UG and[9l Similarly, for A/q = 3, the shape resembles 
the loop turn as seen in Fig. [6l 

In Fig. [3] the curvature-based harmonic- wave descrip- 
tion dHI is compared with the real-space harmonic wave 

y = A' sin((7'x + (/)') (9) 

that was used in [28] to describe worm trails with shal- 
low turns. For low amplitudes of the curvature wave 
A/q < 1 the two relations ([8]) and ([9j) yield nearly equiv- 
alent shapes, but the results significantly differ outside 
this range. An analysis of the skeletonized worm shapes 
depicted in Fig. 3] shows that our harmonic-curvature 
model ([8]) describes nematode body postures accurately 
for all amplitudes. The real-space harmonic description 



([9]) becomes inaccurate already for moderate-amplitude 
undulation (cf. Fig. |4p). 



Piecewise- harmonic-curvature function 

Our simple harmonic-curvature model (|8|) can be ap- 
plied to describe both individual worm postures and 
entire nematode trails (such as the trail depicted in 
Fig. [1]). Complete trails that include gradual and deep 
turns can be represented with high accuracy using a 
piecewise-harmonic-curvature function. We first discuss 
piecewise-harmonic representation for individual nema- 
tode postures, and then we analyze trajectories of crawl- 
ing worms. 

Worm body postures. The shape of C. elegans 
shown in Fig. [5] has a distinctly higher curvature in the 
head section than in the tail section. For such a shape, a 
single-mode sinusoidal-curvature representation is insuf- 
ficient. Figure [5^ shows a line obtained by integrating 
Frenet-Serret equations (|6|) with the harmonic curvature 
([8]) and parameters adjusted to fit the tail part of the 
worm. From the tip of the tail up to the point indicated 
by the filled circle, the line matches the skeletonized ne- 
matode posture very well. Beyond this point a sudden 
deviation occurs, as determined by the local error shown 
in the inset. We conclude that the worm exhibits more 
than one mode along its body length, and there is a sud- 
den change in the curvature of the worm's body at the 
point marked, where the transition occurs. 

To accurately describe the entire worm posture, we use 
the piecewise-harmonic curvature 

|^A2 Sm((725 + 02); Si < S < Su- 

Here 5i is the intermediate point at which a sudden tran- 
sition between harmonic modes occurs. As demonstrated 
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FIG. 7. Piecewise-harmonic-curvature (PHC) representation of a C. elegans trail. (A) Worm trail obtained by 
superposing individual video frames. (B) Independent single-mode harmonic-curvature fits (solid lines) to different regions of 
the skeleton data (open circles). Different segments are shown in different colors, and the filled square marks the beginning 
of the trail. Sudden deviations of the single-mode fits from the skeleton data indicate that the worm switched its mode. The 
insets show the local fit error (pT)) for the second and fourth trail segment. (C) A continuous PHC representation obtained by 
fitting the solution of equations (|6]) and to the trail-skeleton points. The end of each mode is marked with a solid circle. 
The inset shows the local fit error. 



in Fig.[5p, the two- mode curvature function yields an ac- 
curate description of the shape of the whole worm from 
tail to head. 

Two-mode fits of similar quality have been obtained for 
other nematode shapes, such as the loop posture depicted 
in Fig. [6l At the switching point si the worm curvature 
is often discontinuous, but the worm contour obtained 
via integration of the second-order Frenet-Serret differ- 
ential equation (j6j) is continuous and has a continuous 
orientation. 

Worm trails. For C. elegans crawling without slip 
on a solid surface, the instantaneous body postures can 
be treated as segments of the entire worm trail (cf. Fig. 
[1]). Therefore, we assume that the abrupt changes in 
the parameters of the curvature wave that propagates 



throughout the worm body are reflected in the overall 
trail shape. 

A typical example of a trail made by C. elegans crawl- 
ing on an agar surface is depicted in Fig. [TK- The picture 
was obtained by superimposing consecutive frames of a 
video recording of a crawling worm, as described in the 
Methods section. The worm undergoes a small amount 
of slip, which is reflected in occasional trail imperfec- 
tions. These imperfections are also due to the random 
movement that the head of the nematode makes before 
deciding on the path of its motion. However, once the 
head continues to move in a certain direction, the rest of 
the body follows the head on a determined path. 

The trail of a gradually turning worm shown in Fig.[7K 
cannot be represented using a single harmonic-curvature 
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FIG. 8. Parameters of the PHC representation. (A) 

The normahzed amphtude A (blue - sohd), wavevector q 
(green - dashed), and phase (red - dotted) of the curva- 
ture wave along the trail shown in Fig. [71 (B) The evolution 
of curvature along the same trail. 



mode. However, similar to the results for worm postures 
(cf. Fig. [5]), individual segments of the trail are consis- 
tent with the harmonic-curvature representation, as il- 
lustrated in Fig. [7B. In a certain range of s each local 
harmonic-curvature fit follows the track with high accu- 
racy, and then rapidly deviates from the track as shown 
in the insets in Fig.[7B. The errors of fits are less than 1% 
up to the mode transition point and then show a steep 
increase. 

Based on the above observations, we expect that the 
entire trail can be well represented by a continuous line 
corresponding to the piecewise-harmonic curvature 



' A\ sin(g'i5 - 
A2 sin(g'25 - 



>2); 



^ Ai sin((7iS + ^i); 



8^ < S < S\ 
Si < S < S2 

Si < S < Se. 



(11) 



Here 5i, 52, ... Si denote the intermediate points of the 
worm's trail where a sudden transition between modes 
occurs, and sq and Se are the beginning and end points 
of the trail, respectively. The continuous PHC represen- 
tation for the whole trail is depicted in Fig.[7p. The inset 
in Fig. [7p shows the local deviation of our PHC model 
from the skeleton data. The low deviation (less than 1% 
with occasional jumps to about 2%) for all values of the 



arc-length parameter s indicates that our model fits the 
data with high accuracy. 

Figure [8]A. shows the discontinuous variation in the 
curvature- wave amplitude, wavevector and phase along 
the worm trail depicted in Fig. [TK- Jumps of these three 
parameters occur at irregular intervals, and the size of 
the jumps is also irregular. The corresponding evolution 
of the curvature is illustrated in Fig. [Hp. The curvature is 
discontinuous at each mode shift, but the final real-space 
PHC representation, obtained via integration of the sec- 
ond order differential equation (|6]), is continuous. 

Further examples of worm trails and the correspond- 
ing PHC representations are shown in Figs. [9] and [TOl 
Our description captures well both the gradual changes 
of the overall path direction (cf. Figs. [71 and [TO]) and rapid 
changes such as ^l-turns and other deep turns (cf. Fig. [9]). 
An analysis of the variation in the curvature-wave pa- 
rameters indicates that a deep turn may result from an 
increase of the wave amplitude, a decrease of the wave 
vector, or from both such changes occurring at the same 
time. 

The nematode C. elegans changes suddenly its undu- 
lation mode not only during forward motion, but also in 
the backward crawling gait. It has been observed [29| 
that upon reversal of its direction of motion C. elegans 
often retraces part of its track, and then abruptly veers 
from the previous path. This deviation has a very similar 
geometry to the abrupt divergence between the trail of a 
worm crawling forward and a single-mode fit to the trail 
(cf. Fig. [7B). We conclude that both deviations result 
from the same general behavioral pattern, i.e., from a 
sudden change in parameters of the harmonic-curvature 
wave propagating through the body of C. elegans. 



DISCUSSION 

Our piecewise-harmonic-curvature model provides a 
simple analytical description of the crawling motion of C. 
elegans in terms of elementary sinusoidal movements in 
curvature space. We have demonstrated that the worm 
turns by changing the curvature- wave parameters. In 
this section we describe consequences of our findings and 
indicate directions for further research. 



Curvature evolution in the intrinsic worm 
coordinates 

Our observations provide important clues on how the 
neural system of C. elegans controls the nematode move- 
ment to generate propulsion and navigate the envi- 
ronment. To elucidate this aspect, we transform our 
curvature-based description to the intrinsic worm coor- 
dinates with the arc length = s — ST{t) measured from 
the position of the tip of the nematode tail, instead of the 
beginning of the trail. In this representation, the curva- 
ture at position s^ and time t can be expressed in the 
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FIG. 9. PHC representations of trails and the corresponding curvature-mode parameters. Left column: experi- 
mental images of several trails of C. elegans. Middle column: the PHC representations ()lip (solid lines) of the skeletonized trail 
data (open circles). The beginning of the whole trail is marked with a solid square; the end of each curvature mode is marked 
with a solid circle. Right column: The normalized amplitude A (blue - solid), wavevector q (green - dashed), and phase (red 
- dotted) of the curvature modes along the trail. 



form 

K^{s' ,t) = n{s' ^Vst), 0<s<L, (12) 

where is the piecewise-harmonic function (fTTj) . To de- 
scribe the tail position s^it) we used Eq. (|3]), on the 
assumption that the worm moves along the trail with a 
constant velocity Vs. 



The worm's point of view on curvature modes. 

According to our model, the worm motion is entirely de- 
termined by the initial worm configuration at t = 0, 

K^{s', 0) = k{s'), 0<s' <L, (13) 

and the boundary condition 

K.^{L,t) = K.{L^Vst), t>0 (14) 
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FIG. 10. PHC representations and the corresponding curvature-mode parameters for additional trails. 



that describes the head curvature at time t. The curva- 
ture generated in the head segment propagates along the 
worm body according to relation ([12]). 

It follows that the worm (i.e., the nervous system that 
controls its body movements) needs to perform three, 
relatively simple, independent tasks: 

(i) generate a harmonic oscillation of the curvature in 
the head part of the body, corresponding to the 
relation 

n^{L,t) = A^m{ujt^(l)'), (15) 
where uo = qVg is the angular frequency; 

(ii) propagate the curvature wave backwards along the 



nematode, as described by Eq. (p!2]) : 

(iii) produce random jumps of the amplitude A, 
wave vector and phase (j)' of the harmonic time 
variation in the head curvature ([15]). 

Since the curvature is an intrinsically local quantity, all 
three tasks require only local operations. Generation of 
the curvature in the head segment (tasks (i) and (iii)) 
involves local muscle contractions. Propagation of the 
curvature wave (task (ii)) depends on local muscle con- 
tractions and communication between neighboring seg- 
ments. 

The simplicity of our geometrical description strongly 
suggests that there is an underlying simphcity in the way 
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FIG. 11. Body postures of C. elegans swimming in 
water. (A) C-shaped and (B) sine-shaped nematode pos- 
tures. Left column: experimental images. Right column: the 
corresponding skeletonized data and single-mode harmonic- 
curvature fits. 



C. elegans controls its motion. The worm is not aware of 
its entire posture, which would require long-range com- 
munication between spatially separated body segments, 
so the nematode controls its body locally The above 
conjecture is also consistent with results of a recent study 
by Stirman et al. [30|. By targeted illumination, these 
researchers have manipulated the behavior of optogenet- 
ically engineered C. elegans. In one of their experiments 
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FIG. 12. The phase evolution cj) for a sequence of im- 
ages of C. elegans swimming in water. The sequence of 
images from the movie of a swimming C. elegans was fitted 
using harmonic-curvature function (|8]) with the phase cj) vary- 
ing from frame to frame and the amplitude A and wavevector 
q identical for all frames. 



(cf., Supplementary Video 1 of [SO^) neural excitation 
produced by a laser-light impulse applied to the head 
section of a crawhng nematode created a sharp bend in 
the worm body After this initial perturbation, the bend 
propagated backwards along the body of C. elegans^ caus- 
ing the nematode to follow the direction of the head, and 
resulting in a worm crawling along an artificially prede- 
termined path (in the shape of a triangle). This behavior 
supports our conclusion that the propagation of the cur- 
vature wave along the worm body is independent of the 
curvature-generation mechanism. 



Worms moving in different environments 

In this paper we have focused on C. elegans moving 
without slip on a solid substrate. However, based on our 
discussion of the way the worm controls its body move- 
ments, we expect that key elements of our description are 
likely to apply more generally, i.e., even when C. elegans 
undergoes a significant slip, and the nematode trail is not 
so well defined. 

A typical example of such a situation is C. elegans 
swimming in water. During swimming, the nematode 
undergoes a significant translational and rotational drift 
with respect to the surrounding fluid. Thus a superpo- 
sition of subsequent nematode images does not show a 
well-defined trail. To describe the motion of C. elegans 
in its swimming gait, it is more convenient to analyze 
a sequence of individual images (rather than their su- 
perposition) using intrinsic worm coordinates, i.e., the 
framework defined by equation (p!2|) . 

Two images from a video of a swimming C. elegans are 
shown in Fig. [11] together with the harmonic-curvature 
fits to the numerical skeletonization data. The nema- 
tode posture in each frame can be accurately represented 
by contours corresponding to the harmonic curvature (|8|) 
with a single mode of the same amplitude and wavevec- 
tor A/q ^ 0.96 and qL ^ 4.83, but a varying phase (j). 
This suggests that the internal kinematics of a swimming 
and crawling C. elegans are quite similar. The evolution 
of phase ^ as a function of time for a sequence of video 
images is shown in Fig. [121 The gradual drop in (j) with 
constant A and q suggests that a single wave propagates 
backward as the nematode swims. 

A more detailed analysis is needed to determine 
whether C. elegans in the swimming gait uses a simi- 
lar strategy to turn as in the crawling gait. However, 
based on previous reports that the worm uses ^l-turns 
for changing swimming direction [31], we postulate that 
a swimming worm, similar to a crawling worm, abruptly 
varies the harmonic-wave parameters at irregular inter- 
vals, which results in a modification of the swimming 
course. 

We believe that the framework based on harmonic- 
curvature- wave assumptions ([TT]) and ([T2|) will be useful 
in the analysis of worm propulsion and maneuverability 
in water. We also expect that a similar framework can 
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be applied in microstructured environments, such as soil 
or arrays of microfluidic pillars. 

Conclusions and future directions 

Our study has demonstrated a significant geometrical 
simplicity of the crawling gait of C. elegans. We have 
shown that C. elegans propels itself using a simple set 
of elementary movements, and that both shallow turns 
and sharp loop- and 1^-turns correspond to the same ge- 
ometrical family of contours generated by a harmonic 
variation of the curvature along the body of the nema- 
tode. The geometrical simplicity revealed by our analysis 
directly reflects key underlying aspects of the neural con- 
trol of the nematode motion (as explained above). 

Our analytic description of the nematode kinematics 
provides a powerful tool for quantitative investigations of 
nematode behavior. For example, behavioral changes of 
C. elegans in response to chemical signals can be quan- 
tified by the statistics of jumps of the curvature- wave 
parameters A, g, and (j). Since the worm controls these 
parameters directly (in contrast to the angles of the corre- 
sponding turns), further studies based on our description 
are likely to reveal important information on underlying 
mechanisms of chemotaxis. 

We expect that our analysis will facilitate development 
of advanced models of neural control of nematode motion. 
Current models [32] capture gross qualitative features of 
the crawling gait of C. elegans fairly well. However, 
as implied by our results, more comprehensive models 
should include separate neural-network functional units 
that would (i) control curvature-wave generation in the 
head segment of the worm, (ii) propagate the wave along 
the worm body, and (iii) change the wave parameters. 
Introduction of such models would be an important step 
towards back-engineering of the neural control system of 
C. elegans. Such models would also aid development of 
a control system for future synthetic worm-like particles 
that can crawl autonomously. 

MATERIALS AND METHODS 
Experimental Details 

Worm Preparation. C. elegans were cultured on a 
60x 15 mm petri dish containing 4 wt7o agar at 20° C. The 
worms were fed every 4 days with 0.05 ml of E. coli, and 
the petri dish was subsequently wrapped with paraffin 
film to prevent the growth of other bacteria or mold. To 
conduct experiments, a chunk of agar containing worms 
was transferred onto a fresh bacteria-free 4 wt% agar 
plate js^- Images were acquired once all the worms had 
left the chunk and had begun crawling on the fresh plate. 

Image Acquisition. Images of C. elegans crawling 
on agar plates were taken using a Zeiss Stemi 2000- 
C stereo microscope with a Stream View CCD camera 



(640x480 pixels with a 7.4 fim^ pixel size and an 8-bit 
mono sensor), typically acquiring between 10-30 frames 
per second, and were saved as black and white avi movie 
files. The agar plates were illuminated from below by 
mounting the plates on one face of a right angle prism 
with an aluminized hypotenuse (50 mm leg length, 70.7 
mm hypotenuse length). A Chiu Technical Corporation 
Lumina FO-150 fiber optic illuminator with a focusing 
lens and light diffuser were used to illuminate one right- 
angle face of the prism; the light was directed through 
the sample by the hypotenuse face of the prism, and was 
acquired by the stereo microscope. Worms imaged with 
this apparatus appear as a dark feature on a light back- 
ground. 

Image Analysis. Individual worm shapes were ob- 
tained from the avi movie files by converting the movie 
to 8-bit black and white tiff images using ImageJ soft- 
ware (http://rsbweb.nih.gov/ij/). The individual images 
were then analyzed as follows. A grey scale threshold was 
applied to the images to separate the worm from the 
background, small non-worm objects were removed from 
the images based on their size, 'holes' in the thresholded 
worms' bodies were filled and the solid binary images of 
the worms' bodies were transformed into a skeleton run- 
ning down the center of the body by applying a morpho- 
logical thinning operation. These operations were carried 
out using the Matlab Image Processing toolbox routines 
bwareaopen.m^ imfill.m^ and bwmorph. m lespectWely. In- 
dividual data points were collected from the resulting bi- 
nary skeleton image by determining which pixels had a 
value equal to 1, collecting their row and column posi- 
tions, and utilizing these positions as Cartesian (x, y) co- 
ordinates in order to perform curvature fits. In the case of 
Fig. [6l the above described method of automatic skele- 
tonization using the image processing software did not 
provide analyzable skeleton data due to the overlapping 
body segments of the worm. Therefore, in this case, the 
image processing software was used to obtain the skele- 
ton of the worm except in the overlapping region. For 
the overlapping region, the position of the pixels along 
the centroid of the worm were selected and joined with 
the rest of the worm body. 

Worm trajectories, such as those shown in Figs. [71 
[9l and[TOl were constructed by opening an avi movie of 
crawling worms in IMAGE J, and projecting the darkest 
pixels from the time-series of images into a single image 
(using ImageJ 's Z project tool). This resulting image 
was then processed in the same manner as individual 
worm images in order to obtain (x, y) data points along 
the worms' trajectories. 



Modeling details 

Solution of differential equations. The set of cou- 
pled differential equations (|6]) was solved using the Mat- 
lab's fourth-order Runge-Kutta ODE solver ode^S.m. 
Since Eq. (j6|) is a set of second order differential equa- 
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tions, it requires two initial conditions — the starting 
point and the initial direction. These two conditions were 
obtained by the fitting procedure described below. For 
the continuous representation with multiple modes along 
the trajectory, each mode was determined using the loca- 
tion and direction of the last point of the previous mode 
as the initial conditions. With sudden jumps in param- 
eters for the harmonic-curvature function, the curvature 
is discontinuous between each mode as shown in Fig. [Hp. 
However, the corresponding PHC representation of the 
trajectory in real (Cartesian) space is a continuous curve 
with a continuous slope, because it is a solution of a set 
of second-order differential equations. 

Measure of error. The deviation of the model from 
the experimental data was measured using the error func- 
tion 

e, = L-'(^l±d^\ (16) 

where n is the number of data points, and di is the dis- 
tance of the curve ^ from the data point i. The quantity 
(p!6|) was minimized to obtain the model parameters for 
accurate description of the worm body postures and tra- 
jectories. In all cases except one we used p = 2. For 
the trail presented in Fig. [9K, we used p = 12 to avoid 
large local deviations between the model and the data. 
The error (p!6|) was minimized using the built-in Matlab 
optimizer fminsearch.m. 
The local errors 

e{sj) = L-'(^^ > (17) 

shown in the insets of Figs. [5] and [7] were estimated by 
calculating the average distance of the curve from the 
data along a moving window of five data points. 

Model parameters Ai, qi, and The model pa- 
rameters for a given trajectory were evaluated in two 
stages. In the first stage, individual segments of the tra- 
jectory were fitted independently using the single- mode 
representation ([8]) to find the positions of the mode shifts 
and to determine the initial guesses for the parameters 
Ai^ Qi^ and In the second stage, a continuous multi- 
mode representation of the trajectory was obtained using 
the initial guesses from the first step. 

Single-mode independent piecewise fits. The fitting 
procedure was started from one end of the trail with ap- 
proximately 25 consecutive skeleton points. The model 
parameters for this section were obtained by minimizing 
Eq. (p!6]) . The quality of the fit was estimated by cal- 
culating the local deviation between the model and the 
experimental data using Eq. (pT|) . If the local error of a 
particular segment showed no increase, the length of the 
section was doubled and the above steps were repeated. 
If the local error showed a steady increase beyond a cer- 
tain point, as seen in Fig. [7p, the fit was truncated at 



that point and the point of truncation served as the lo- 
cation of a mode shift. The above steps of independent 
harmonic fits were repeated until the end of the trajec- 
tory was reached. 

Continuous Fit. The purpose of this stage was to de- 
termine the parameters A^, Qi^ and for the continu- 
ous representation of a trajectory according to the PHC 
model (pT]) . The continuous representation was obtained 
by consecutively fitting larger and larger portions of the 
trajectory to avoid difficulties associated with local min- 
ima of the accuracy measure (pT|) . The fitting procedure 
was started from the first segment and involved addi- 
tion of subsequent segments in two steps. In the first 
step, the curvature parameters of segment k were refit- 
ted with the constraints of maintaining the continuity of 
the curve and its slope between the current and previous 
segments. In the second step, segments 1 through k were 
readjusted allowing variations of their curvature param- 
eters and locations of mode shifts. Typically, we allow 
only variations of parameters of segments A:-2, A;-l, and k 
in this step. The above two steps were repeated until the 
end of the trajectory was reached. 

In several cases (Figs. [71 [9p, and [TOK). the above 
procedure did not provide an accurate match between 
the model and the trajectory. This difficulty was asso- 
ciated with defects of the trail due to the local slip of 
the worm. For these cases, the continuous PHC repre- 
sentation of the trajectories was obtained by first fitting 
the two sections (before and after the region with imper- 
fections) independently, and then joining them smoothly 
using an additional segment corresponding to the region 
of the imperfect data. 

In all cases, we achieved the accuracy e {sj) ^0.02 and 
the typical segment length was 0.7 L. The fits performed 
in the forward and reverse directions yielded similar pa- 
rameter values, which indicates that our model is robust. 
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